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An experimental study was conducted to determine the proximity aeroheating character- 
istics for a two-stage-to-orbit concept in close proximity in the NASA Langley 20-Inch Mach 
6 Air Tunnel. A new hybrid discrete thin-film resistance gauge technique was evaluated in 
this study and used to measure experimental interference heating levels between the booster 
and the orbiter at a constant freestream Reynolds number of 8.25x1 0 6 /m and a variety of sep- 
aration and axial offset distances. It was found that, as the orbiter separates from the booster 
and the booster falls away, the windward centerline heating increased on the orbiter by as 
much as 13-times over the baseline, single model heating distribution, and on the booster by 
as much as 6-times. The aeroheating database developed can be used for computational fluid 
dynamic code validation. 


Nomenclature 

h = heat transfer coefficient, h=q/(H aw -H w ), (kg/m 2 /s) 

H = enthalpy (J/kg) 

L = model reference length (m) 

p = pressure (Pa) 

q = surface heat transfer rate (W/m 2 ) 

R n = model reference nose radius (m) 

Re = unit Reynolds number (1/m) 

T = temperature (K) 

U = velocity magnitude (m/s) 

x = axial distance from nose of model (m) 

y = spanwise distance from centerline of model (m) 

a = angle-of-attack (deg) 

p = density (kg/m 3 ) 

Subscripts 

oo m freestream static conditions 

aw = adiabatic wall conditions 

FR = conditions from Fay-Riddell calculation for a hemisphere 
tl = reservoir conditions 

t2 = stagnation conditions behind a normal shock 
w - wall conditions 
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I. Introduction 


N ASA’s Integrated Space Transportation Plan (ISTP) provided a comprehensive, long-term strategy to meet 
future launch vehicle and technology needs combined with reduced cost and improved safety and reliability. 
Given the technology levels in the areas of propulsion, structures, and materials, the technical risk for Single-Stage- 
To-Orbit (SSTO) was high, and advanced Reusable Launch Vehicle (RLV) design attention shifted back to Two- 
Stage-To-Orbit (TSTO) concepts. While TSTO concepts which were considered ranged greatly with respect to 
orbiter/booster shape, propulsion system characteristics, and resulting nominal flight trajectories, all TSTO configu- 
rations presented a critical need to enhance the aerodynamic and aerothermodynamic knowledge base in two-body 
separation characteristics within the sensible atmosphere to satisfy both nominal and abort ascent separation scenar- 
ios. The present work adds experimental aeroheating data at Mach 6 to the existing database of aerodynamic data pre- 
sented in Ref. 1. This database of an orbiter and booster of comparable size in close proximity can be used as a 
computational fluid dynamics validation tool. A new hybrid thin-film resistance gauge technique was evaluated dur- 
ing the present study to measure aeroheating levels on vehicle surfaces where optical access is limited and the tech- 
nique of phosphor thermography could not be used. While this effort was initiated to satisfy specific requirements in 
the former Next Generation Launch Technology (NGLT) program 2 , the tools and methodologies developed were 
generic in nature and may be applicable to a broad range of launch vehicles/architectures in NASA’s new Space 
Exploration Initiative. 


II. Configuration Description 


The same generic two-stage configuration that was used 
in the aerodynamic study in Reference 1 was also used in the 
present study. The Langley Glide-Back Booster (LGBB), in a 
belly-to-belly, bimese configuration is shown in Fig. 1. The 
LGBB is a generic wing-body configuration that was devel- 
oped by the Langley Research Center Vehicle Analysis 
Branch for in-house multi-stage launch vehicle system-analy- 
sis studies. 

Figure 2 shows a detailed planform schematic of the 
LGBB. The full-scale configuration has an 2.09-m diameter 
cylindrical fuselage with a drooped, 0.28 m radius nose. The 
length from the tip of the nose to the base of the body (i.e. the 
engine exit plane) is 19.05 m. The wing/strake is a NACA 
2300 series airfoil at a 6-degree dihedral angle; the strake is 
swept 8 1.1 -degrees and the wing is swept 45.3-degrees. The 
wing’s forward swept trailing edge is a full span elevon pair 
separated by a center body flap. The configuration has a cen- 
ter vertical tail with a leading edge sweep angle of 45- 
degrees. The tail’s airfoil section, including a trailing edge 
rudder, is a NACA 0013 airfoil and its span is 2.80 m. The 
baseline LGBB vehicle includes a swept canard, NACA 
64A012 airfoil, located on the forward part of the fuselage. 
The vertical tail and canards, however, were not included in 
the geometry tested in the present study to simplify the appli- 
cation of the thin-film gauges (discussed below). 

III. Experimental Methods 
A. Test Facility 

The 20-Inch Mach 6 Air Tunnel is a blowdown facility in 
which heated, dried and filtered air is used as the test gas. The 
tunnel has a two-dimensional, contoured nozzle which opens 
into a 0.5207 m by 0.508 m test section. The tunnel is 



Figure 1. LGBB bimese arrangement 



Figure 2. Dimensions of the LGBB vehicle. 
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equipped with a bottom-mounted injection system that can transfer a model from the sheltered model box to the tun- 
nel center-line in less than 0.5 seconds. Run times of up to 15 minutes are possible in this facility, although for the 
current aeroheating study run times of only a few seconds were required. The nominal reservoir conditions of this 
facility are stagnation pressures of 206.8 kPa to 3447.4 kPa with stagnation temperatures of 422.2 K to 555.5 K, in 
which the freestream and post-shock flow behaves as a perfect gas (y ~ 1.4) with Mach numbers between 5.8 and 6.1 
and Reynolds numbers of 1.64xl0 6 /m to 23.95x1 0 6 /m. The nominal flow conditions for this study are presented in 

Table 1. A more detailed description of this facility is presented in Ref. 3. 

Table 1. Nominal flow conditions for the 20-Inch Mach 6 Air Tunnel for present study. 



Re<x, (1/m) M, x p t (N/m 2 ) 

T,( K) 

Poo (kg/m 3 ) Uao (m/s) H, (J/kg) h FR 

8.2789E+06 5.9511 1.0533E+06 

503.15 

3.9408E-02 942.23 5.0612E+05 0.9413 



B. Test Model Description 

Cast ceramic test models were used in the present study due to the speed and cost effectiveness at which they can 
be produced. In order to manufacture the ceramic test models, rapid-prototype, stereolithographic resin models were 
first fabricated based on surface geometry definitions in electronic data files. Wax molds of the resin models were 
made, and then a patented 4 silica ceramic slip casting technique was used to form ceramic shells of the models. The 
shells were then back-filled with a hydraulically setting magnesia ceramic for strength and support. The cast ceramic 
aeroheating models were 0.0 121 -scale representations of the LGBB (models were 0.23 m long). 

C. Test Techniques 

Although the primary test technique of the present study was a hybrid thin-film resistance gage technique, lim- 
ited measurements were made using the technique of global phosphor thermography. These techniques are described 
below: 

1. Phosphor Thermography Technique 

For two of the cast ceramic models, global surface heating distributions were measured using the global optical 
method of two-color, relative-intensity, phosphor thermography 5 " 8 . The ceramic wind tunnel models were coated 
with a phosphor compound that fluoresces in two separate regions (green and red) of the visible light spectrum. Dur- 
ing a wind tunnel run, the phosphor-coated model is illuminated by ultraviolet (UV) light sources, and the resulting 
fluorescent intensity of the model is recorded and digitized through a color CCD (charge coupled device) camera. 
The fluorescent intensity is dependent on both the intensity of the incident UV light and the local model surface tem- 
perature. The incident UV intensity dependence is removed by taking the ratio of the red to green intensity images, 
from which surface temperature distributions can be determined through prior calibrations. Images are acquired 
before the wind tunnel run and after injection of the model to the tunnel centerline during a run. Global heat transfer 
distributions are then computed from these temperature data using one-dimensional, semi-infinite solid, constant 
heat- transfer coefficient conduction theory 8 . 

2. Thin-Film Resistance Gauge Technique 

The global phosphor thermography technique could not be used to measure heating rates between the two vehi- 
cles because of line of sight restrictions. Therefore, a thin-film instrumented model was fabricated in a similar man- 
ner to that of Ref. 9. As stated above, cast ceramic models, usually used in phosphor thermographic aeroheating 
studies, were used in the present study. A thin polyimide film was then used as a transfer medium to place the gauges 
down on the vehicle surface. The thin-film gauges were etched onto the flat polyimide film, and the film was later 
bonded to the vehicle surface. The polyimide film that was chosen was a 5. 08x1 0' 5 m thick type S Upilex® film. 
(Upilex is a registered trademark of Ube Industries, Ltd.) 

The thin-film gauge layout is presented in Fig. 3. A long row of 44 gauges, spaced 5.08 mm apart, was placed 
along the vehicle centerline. Gauges were also placed symmetrically about the centerline on the upper and lower sur- 
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faces of the wings. On each surface of the wings, a total of 13 gages were placed in 3 rows. A total of 96 gages were 
applied to the thin- film model. 




Windward Surface Gauges. 


Leeward Surface Gauges. 


Figure 3. Thin-film gauge layout. 

To facilitate the application of the gauges onto the ceramic model, the film panel was split into five separate pan- 
els (centerline and each surface of the wings). An adhesive was then used to bond the panels to the ceramic substrate 
one at a time. Each panel was positioned on the model, which was then placed in a plastic vacuum bag to mold the 
film to the model surface. As the air was suctioned from the bag, the panel was monitored to make sure it was located 
in the correct position and orientation. After this portion of the panel application was complete, the model with the 
newly applied panel was taken out of the vacuum bag and the panel was smoothed further to take out any adhesive or 
air bubbles under the panel. Any excess adhesive was removed. 

The application of the film panels turned out to be more of an art than a science. The process of applying all five 
panels onto the model resulting in a satisfactory surface fidelity (no gaps at the split lines or waves in the panels) took 
approximately eight months since unsatisfactory panel applications had to be removed from the ceramic surface, 
which was carefully cleaned, and reapplied until the panel bond was acceptable. 


D. Data Reduction 

1. Phosphor Thermography Technique 

One-dimensional, semi-infinite solid heat conduction theory 8 was used to compute surface heating distributions 
from the global surface temperature data acquired through the technique of two-color, relative-intensity, phosphor 
thermography. A constant heat-transfer coefficient is assumed in this theory, and empirical corrections 8 are made to 
account for changes in model substrate thermal properties with temperature. Phosphor images were acquired shortly 
after injection of the model to the tunnel centerline, which requires less than one second. 

Data were extracted in the IHEAT 7 software package axially along the centerline. Results are presented herein in 
terms of a non-dimensional heat transfer coefficient ratio, h/h FR , where h FR is the theoretical stagnation point heating 
to a 3.44xl0" 3 m radius sphere (the nose radius of the test models) computed using the method of Fay-Riddell 10 for a 
wall temperature of 300 K. 

2. Thin-Film Resistance Gauge Technique 

A calibration was performed on the thin-film model (gauges already applied to ceramic model) over a tempera- 
ture range of 291-395 K. The model was placed in an oven and allowed to come to a steady temperature in approxi- 
mately 10 degree increments. The millivolt output of each gauge was recorded and then plotted vs. the temperature to 
back out a gauge calibration curve fit. Ordinarily, a second-order curve fit is used to convert the millivolt output of 
the gauges to temperature. However, in the temperature range tested in the present study, a linear fit was all that was 
necessary due to the linearity of the millivolt- temperature response. 

Data reduction was performed using the code 1DHEAT 11 . The finite volume implementation of 1DHEAT takes 
direct account of the variable substrate thermal properties and removes the restriction of a semi-infinite substrate 
(thus allowing for a multiple layer analysis). This was important because of the dual layer of Upilex and ceramic on 
this model. Thus, heat transfer was computed using the finite volume technique, with 50 points in the Upilex layer 
(5.08xl0 -5 m thick) and 150 points in the ceramic layer (6.35xl0“ 3 m thick). The current analysis has been conducted 
with the assumption that the approximately 25.4-pm bond-line was negligible. Reliable thermal property information 
for the high-temperature epoxy used was not available. 
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E. Error Analysis 

1. Phosphor Thermography Technique 

The estimated experimental uncertainty in heat transfer of the thermographic phosphor system is a function of 
model injection time, substrate thermal properties, initial wall temperature, run-time wall temperature, and fluores- 
cence intensity. For higher model surface temperatures (greater than 400 K), such as those on most of the windward 
surface, the total (bias and precision) uncertainty 8 is approximately ±8% to ±10%, while for lower temperatures (less 
than 325 K), the total uncertainty is approximately ±15% to ±20%. Additional measurement uncertainty can be intro- 
duced due to heat conduction in high-gradient regions such as the leading edge of the LGBB. 

Uncertainties can be introduced when extracting line cut data from an image due to perspective distortion of the 
image, lack of pixel resolution in high-gradient regions, and lack of precision in locating fiducial marks. These errors 
are estimated to be less than ±5% on relatively flat surfaces and up to ±10% on highly curved surfaces. 

A square root of the sum-of-the-squares estimate for the total uncertainty based on the above factors gives a 
worst-case experimental uncertainty range of ±13% on flat areas experiencing high heating to ±25% on curved areas 
with little heating. Since the test model in the present study is at a = 0-deg. for all runs (except for two test technique 
comparison runs), the total uncertainty of the phosphor thermography technique for the present study is on the order 
of ±15%. 

2. Thin-Film Resistance Gauge Technique 

Primary contributions to the uncertainty in the heat-transfer results for the thin-film resistance sensors include 
the uncertainties associated with the thermal properties of the substrate material. This includes both the neglected 
bond layer and the porous ceramic sub-layer which absorbed some of the adhesive. For traditional thin-film studies 3, 
11_13 , the data is believed to be accurate to within 8-10% in regions where one-dimensional heat conduction can be 
assumed. The uncertainty levels for the thin-film resistance gage model of the present study should be higher than for 
the traditional techniques, as the thermal properties of the bond line and ceramic/adhesive sub-layer are not accu- 
rately known. Assuming the bond layer introduced a 1% error and the thermal properties of the ceramic sub-layer did 
not change by any more than 15%, a square root of the sum-of-the-squares estimate for the total uncertainty gives 
±20%. 

F. Test Matrix and Tunnel Conditions 

The data were collected on the LGBB models at angles-of- attack of 0-, 20-, and 3 0-deg with a singe model for 
comparison between the two test techniques and at a = 0-deg for the bimese configuration. In the bimese configura- 
tion, a phosphor model was tested opposite to the thin-film model. Two geometric separation parameters were tested 
in the present study (see Fig. 4): Az was the distance between the lower surfaces of the vehicles; Ax was the axial dis- 
tance between the noses of the vehicles. All geometric separations are given with respect to the thin-film model (i.e. 
if the thin-film model is 5 -in. behind the phosphor model, Ax = -5 -in.). Therefore, if the value of Ax is negative, the 
centerline data presented is for the booster. In the present study, zlz varied from 0.5-in. to 3-in. and Ax varied from -5- 
in. to ±6-in (both positive and negative values of Ax were tested to obtain both booster and orbiter heating data). A 
constant Reynolds number of 8.28xl0 6 /m was run to be consistent with the Reynolds number tested in Ref. 1 in the 
20-Inch Mach 6 Air Tunnel. A total of 90 runs were performed in a two week tunnel entry. 



Thin Film Model 


Figure 4. Model separation parameters. 
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IV. Experimental Results 


A. Test Technique Comparison 

A comparison between the test techniques of phosphor 
thermography and the hybrid thin-film resistance gage 
technique is included because the thin-film technique used 
is not a mature technology in comparison with standard 
thin-film techniques. This technique was chosen because 
the models could be fabricated quickly and were inexpen- 
sive. The comparison is presented in Fig. 5 for angles-of- 
attack of 0-, 20-, and 30-degrees. The 20- and 30-deg. 
angles-of-attack are included to extend the comparison of 
the two test techniques beyond the low heating levels expe- 
rienced at a = 0-deg to levels comparable to those experi- 
enced due to shock interactions in the present study. The 
thin-film data is approximately 40% lower than the phos- 
phor thermography data. This is believed to be partially 
due to the uncertainty in the substrate thermal properties, 
but was not fully understood as of the writing of this paper. 
This discrepancy is still being investigated. Because of the 
uncertainty in the heating rates, centerline heating mea- 
surements will be presented herein as a ratio of the heating 
rates measured during a run to the heating rates measured 
on a baseline, single model run. Therefore, if the heating 
ratio is equal to 1.0, the heating rate measured at a specific 
gage was equal to the heating rate at the same gauge from 
the baseline run. 
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B. Windward Heating 

Centerline heating augmentations for the bimese con- 
figuration are presented in Figs. 6-17, along with corre- 
sponding schlieren images for the closest separation value 
at each axial offset. This series of figures begins with the 
thin-film model 6-in. in front of the phosphor model (data 
from orbiter) and decreases in increments of 1-in. to the 
thin-film model 5-in. behind the phosphor model (data 
from booster). In each of these plots, the model separation 
(Az) begins at 0.5 -in. and increases in 0.5 -in. increments up 
to 3.0-in. Each plot includes a sketch of the model orienta- 
tion. The thin-film model is labeled as either the orbiter or 
booster in each figure. The legend for this series of figures 
can be seen in Table 2. In general, as the vehicles separated 
(Az increased), the heating augmentation decreased. 

Data and a schlieren image for Ax = 0-in. are presented 
in Fig. 12. A reflected shock train can be seen between the 
two vehicles resulting from the bow shocks. The initial 
bow shock impingement from the opposite vehicle (for the 
closest separation) was located at approximately x/L = 

0.30. The first reflected shock impingement is located at 

approximately x/L = 0.50. Finally, a third shock impingement can be seen at approximately x/L = 0.70. These 
impingement locations correlate well with the schlieren image. 

When considering the orbiter, the initial shock interaction from the booster moved aft with increasing Ax and the 
magnitude of the resulting surface heating increased. The interaction of the bow shocks moved towards the nose of 
the booster, which changed the angle at which the booster bow shock impinged on the orbiter, thus changing the 
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Figure 5. Comparison of test techniques at a = 0- 
20-, and 30-deg. 
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nature of the shock, boundary layer interaction. As the booster bow shock moves aft, the separation region in front of 
the shock impingement can be seen to increase in size (area of decreased heating immediately before increase in heat- 
ing). The largest deviations from the baseline, single model heating distribution occurred at the greater axial separa- 
tions. At Ax = +6-in., heating rates increased by as much as 13-times toward the aft end of the model due to the 
interaction with the bow shock from the booster. It also appears that at the greater axial offsets, the bow shock from 
the booster is also interacting with the shock off of the chine and wing of the orbiter, which could further increase 
heating levels. 

When considering the booster, the initial bow shock impingement moved forward as Ax decreased and the mag- 
nitudes of the impingements decreased. The largest deviation from the baseline heating distribution experienced was 
when the booster and orbiter were even with each other (Ax = 0-in.). In this configuration, both vehicles experienced 
heating augmentations of approximately 6-times what the baseline configuration experienced. As the axial separation 
between the booster and orbiter increased, the heating augmentation decreased over the centerline of the booster. 

Table 2. Centerline augmentation legend. 

Az = 0.5-in. 

Az = 1.0-in. 

Az = 1.5-in. 

Az = 2.0-in. 

Az = 2.5-in. 

Az = 3.0-in. 
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Figure 6. Windward centerline heating augmentation on the orbiter at Ax - +6-in. 
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Figure 7. Windward centerline heating augmentation on the orbiter at Ax = +5-in. 
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Figure 8. Windward centerline heating augmentation on the orbiter at Ax = +4-in. 



0.00 0.20 


0.40 „ 0.60 

x/L 



Figure 9. Windward centerline heating augmentation on the orbiter at Ax = +3-in. 
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Figure 10. Windward centerline heating augmentation on the orbiter at Ax = +2-in. 
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Figure 13. Windward centerline heating augmentation on the booster at Ax = -1-in. 
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Figure 15. Windward centerline heating augmentation on the booster at Ax = -3-in. 
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Figure 16. Windward centerline heating augmentation on the booster at Ax - -4-in. 
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Figure 17. Windward centerline heating augmentation on the booster at Ax - -5-in. 

C. Leeward Heating. 

Global phosphor thermography images of the leeward surface from the Ax = +6-in. bimese configuration are pre- 
sented in Fig. 18 with the corresponding schlieren images in Fig. 19. At this location, the phosphor model was 6-in. 
behind the thin-film model (phosphor model was the booster), so the phosphor model experienced the greatest range 
of possible shock interactions as the thin-film model separated away from it. Most of the leeward surface saw no 
change in heating levels as the booster was moved away from the orbiter. It is believed that at the lower separation 
distances, the orbiter bow shock layer feeds the leeside flow of the booster (Fig. 19. a). At separations of 1.5 and 2.0- 
in., a slight increase in heating was observed along the upper surface of the body near the aft extending from the 
beginning of the wing chine to the aft of the vehicle. It is believed that this increase is due to the interaction of the 
wing and bow shock of the thin-film model, respectively, and the orbiter shock layer no longer feeds the upper sur- 
face flowfield. 




Figure 18. Leeward global aeroheating on the booster at Ax = +6-in. 
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Figure 19. Schlieren images at Ax = +6-in. 

V. Conclusion 

An experimental aeroheating study was conducted in the NASA Langley 20-Inch Mach 6 Air Tunnel on a 
generic bimese two-stage-to-orbit configuration. This test was conducted in part to study the effects of shock interac- 
tions between a booster and an orbiter of comparable size. It was found that, as the orbiter moved ahead and away 
from the booster, heating levels along the centerline of the orbiter were augmented by as much as 13 -times the value 
experienced on baseline, single model. The booster experienced heating augmentations of as much as 6-times the 
baseline, single model configuration. The leeward surface of the vehicles did not experience any significant changes 
in heating levels as the vehicles separated. The database developed can now be used for computational fluid dynamic 
code validation. The present study also evaluated a new hybrid, multi-layer, thin- film gauge technique. It was found 
that the complexities in gauge application and uncertainties in substrate thermal properties outweighed the advan- 
tages of quick model fabrication and lower cost as compared to traditional thin-film techniques. 
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